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Abstract. We compare quiescent and flare X-ray spectra of the RS CVn binary a Gem obtained with the Chandra and XMM- 
Newton grating spectrometers. We find that in addition to an overall 25% flux increase, which can be ascribed to variations in 
the system's quiescence activity over the 15 months that passed between the observations, there is a hot plasma component of 
kT e i 3 keV that arises with the flare. The hot component is manifested primarily by emission from high charge states of Fe and 
by a vast continuum. The cooler (kT e S 2 keV) plasma remains undisturbed during the flare. We find no significant variations 
in the relative abundances during the flare except for a slight decrease (< 30%) of O and Ne. 
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1. Introduction 

The interplay between steady coronal emission and coronal 
flares has been a subject of ongoing research for many years. 
X-ray line resolved spectra available with Chandra and XMM- 
Newton allow now for unprecedented plasma diagnostics, one 
of which is elemental abundance measurements. 

Measurements of relevant changes could offer important 
diagnostics for the heating and the plasma transport pro- 
cess in stellar coronae. It is unclear how exactly large flares 
affect the abundances in stellar coronae. Although indica- 
tions for changes in the metallicity of a flaring corona were 
reported early on fro m low-resolution devices (see discus- 
sion in iGiidel [2004), more reliable analysis had to wait 
for the advent of high-sensitivity and medium-to-high reso- 
lution spectrometers. Indications for increasing metallicities 
were found from ASCA and BeppoSAX medium-resolution 
observations (e.g. iMewe et all 1 19971 iTsuboi et all 1 19981 
Favata & Schmitt , 1999). Moreover, when analysing elemen- 
tal abundances of individual elements, selective, significant en- 
hancements of low fi rst ionization pote n tial (FIP) el e ments 
have been reported (Gud el et all 1 19991: lOsten et all [2000, 
l2002l) . In contrast, in several other cases, the composition of 
the plasma appeared to remain uncha nged during flares (e.g., 
IMaggio et aukoodlFranciosini et alll200lh . 

A crucial reconsideration of the situation came with the ad- 
vent of high-resolution spectroscopy with XMM-Newton and 
Chan dra, but the findings so far still lack systematic trends. 
While Audard 'et all J200ll) found significant enhancements of 
low-FIP elements in a flare on HR 1099, two flares reported 
by Osten et al. (2003; for cr 2 CrB) and Gudel et al. (2004; for 



Proxima Centauri) showed an increase of the abundances of 
several elements, but no selective FlP-dependence was found. 

RS CVn binary systems are bright X-ray and EUV 
sources owing to rapid rotation, generating a magnetic dy- 
namo. As such, they have been studied extensively in both 
bands JAudard et all 120031: Isanz-Forcada etail I2002L respec- 
tively, and references therein). The RS CVn cr Geminorum 
(HD62044, HR 2973, HIP 37629) is particularly bright and 
well observed at all waveleng ths. For an RS CVn, it h as a rather 
long period of 19.6045 days (Duem mler et allll997h . The pri- 
mary star is a Kl III type, red giant. Little is known of the 
secondary as it has not been detected at any wavelength, but 
restrictions to its mass and the low luminosity suggest that it is 
most likely a late-type main sequence star of under one solar 
mass iDuemmle r et all f 1997b . 

The cr Gem system is a luminous X-ray (logLx ~ 31.0 + 
0.2 erg s'^ lYi et'all 1 19971 corrected for a distance of 37.5 pc) 
and radio (logLg at 15.40 erg s Hz -1 , at 6 cm wavelength; 
iDrake et'all Il989t> source. While most observations found it 
to be a relatively steady emitter, a very large flare has been 
detec ted in December 1998 with EUVE JSanz-Forcada et all 
l2002h . Another flare was detected in April 2001 both in the X- 
ray and in the radio in whic h a Neupert effect: ^ t Lx oc LR a dio 
(Neupert, 1968) was found JGudel etall E002). The common 
explanation for the Neupert effect is that the gyrosynchrotron 
radio emission during the flare is due to fast particles spiral- 
ing down the flaring loop. When this particle population hits 
the footprint of the loop they heat the choromspheric material, 
which is subsequently driven into the corona where it emits 
X-rays. Thus, the radio emission depends on the instantaneous 
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Fig. 1. Light Curves for cr Gem in time bins of 500 s. (a) April 
2001 RGS (1 and 2 combined) observation (1st and 2nd or- 
ders), (b) December 1999 LETG/ACIS observation (all orders). 



number of particles running through the loop, while the slowly 
varying X-ray emission from the surrounding chromosphere 
represents the total fast-particle energy converted into heat. 

The detection of the Neupert effect involved only the light 
curve extracted from the EPIC-pn instrument on board XMM- 
Newton. In this paper we investigate the X-ray spectrum emit- 
ted during the 2001 flare and compare it to a quiescent spectrum 
obtained in December 1999. Our main goal is to understand the 
effects of the flare on the thermal and chemical structure of the 
X-ray plasma and to examine whether they support the chro- 
mospheric evaporation scenario. 

2. Observations and data analysis 

2.1. Data Reduction and Light Curves 

The target cr Gem was observed by XMM-Newton in April 
2001 for a total exposure time of 54 ks (observation start: 
2001-04-06 16:24:29, end: 2001-04-07 07:54:40, UTC). The 
data were reduced using the Standard Analysis System (SAS) 
version 6.0.0. In this analysis we use the Reflection Grating 
Spectrometers (RGS) in the 1st order of diffraction, which 
gives reliable data from 6 to 38 A. Line fluxes of Fe 24+ and 
Fe 25+ were extracted from the EPIC-pn data with the use of 
the XSPEC software package ( lArnaudlfl99i) . Background is 
subtracted using off-source CCD regions. Chandra observed 
the target on December 1999 (observation start: 1999-12-27 
18:38:39, end: 1999-12-28T23:07:07, UTC) for a duration of 
100 ks with the Low Energy Transmission Grating (LETG) + 
Advanced CCD Imaging Spectrometer (ACIS) configuration in 
Continuous Clocking (CC) mode. The data were reduced using 
the CIAO package version 3.0.2. 

The long duration fl are, w i th an o bserved peak rise of 20%, 
observed also by Giid el et all ll2002l) . is seen in the RGS light 
curve presented in Figure The rise in flux is roughly uni- 
form in the entire RGS wavelength band. The flare light curve 
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Fig. 2. cr Gem spectra, (a) Flaring ZMM-RGS spectrum, April 
2001, average of the two RGS instruments in 1st order, (b) 
Quiescent state Chandra-LETG spectrum, December 1999, av- 
erage of the two 1st orders. 

is contrasted with the flat light curve obtained from the LETG 
observation shown below in the same figure. 

The use of CC mode is difficult for spectral extraction due 
to order mixing. Above 20 A, the orders become smeared and 
mixed and do not allow for reliable order-sorting. Part of the 
smearing could be due to cosmic ray afterglow, but since there 
is no background assessment for the CC mode, it cannot be sub- 
tracted in a systematic way. Therefore, only data in the range 
of 1.8 to 20 A are used in the present work. One exception is 
the N 6+ line at 24.78 A which can be used. Although the con- 
tinuum level in this region is unreliable, the flux in the line is 
not affected by mixing with the 2nd order, as no major peaks 
occur at 12.4 A. Since the target is very bright, background can 
be neglected. Also, CIAO's default extraction region is too nar- 
row for CC mode, missing out many source photons, especially 
in the -1 order. This is due to CC-mode gain not being well cal- 
ibrated. In order to extract all of the observed source photons, a 
wider energy region had to be selected. In the tg_resolve_events 
task, we used: OSORT_HI=0.34 OSORT_LO=0.20 for m * A > 
7.5 A and OSORT_HI=0.08 OSORT_LO=0.12 for m* A < 7.5 A 
For the spectral analysis, only 1st order data were used. The 
spectra obtained from the quiescent and flare observations are 
shown in Figure |2 

2.2. Direct comparison of flare and quiescence 
spectra 

The most obvious difference between the two spectra in 
Figure |5] is the vast continuum emission in the flare. While in 
quiescence the continuum starts to drop below 12 A , during the 
flare the continuum is still rising at 6 A , which is the limit of 
the RGS. Using the data from EPIC-pn and assuming the con- 
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Fig. 3. Line intesities of the flaring and quiescent states compared. The quiescent lines are multiplied by a factor of 1 .25 and an 
ad-hoc continuum is added to match the flaring continuum. The plotted residual is the quiescent curve subtracted from the flaring 
spectrum, with corresponding error bars. 



tinuum is thermal bremsstrahlung, we can determine the dom- 
inant plasma component contributing to the flare continuum to 
be at kT * 5 keV. 

Given the two high quality spectra, we can compare the line 
fluxes directly and in a model-independent way. In order to do 
that, we need to account for the different line spread functions 
(LSF) of the instruments. LETG has approximately a Gaussian 
LSF with FWHM « 0.05A, while RGS has a slightly wider 
LSF of FWHM « 0.07 A, but a more complex shape with a 
narrow peak. For the comparison, we fold the fluxed LETG 
spectrum through the RGS response. For the details of this pro- 
cess see iNordon et alJ l2QQ$i . In the following, we refer to the 
LETG spectrum, processed through the RGS response matrix, 
as the quiescent spectrum. 

In Figure [3] we compare the flare and quiescent spectral 
lines directly. We multiply the quiescent (LETG) spectrum by 
a factor of 1.25 to match the bright, isolated 19 A 7+ Ly<? line 
flux in the flare (RGS). Subsequently, a smooth phenomenolog- 
ical continuum is added to the rescaled quiescent spectrum to 
match the excess flare continuum. This enables us to compare 
only the line intensities. The flare and processed quiescent (i.e., 
rescaled + added continuum) spectra are plotted in Figure [3] 
Also plotted is the residual of the two with the corresponding 
error bars. 

2.3. Emission measure distribution modeling 

The continuum is almost featureless for a distrubution of 
plasma temperatures and in addition, the continuum from the 
LETGS instrument is unreliable due to the CC mode. Therefore 
line fluxes alone are used to derive the emission measure dis- 
tribution (EMD). Table[2lists the measurable lines in the flare 



and in quiescence. The observed line flux Ft of ion q due to 
the atomic transition j — > i can be expressed by means of the 
element abundance with respect to Hydrogen A z , the distance 
to the object d, the line power P q - t and the ion fractional abun- 
dance f g as: 

A r°° 

F lf=4^2 J P C '/T)MT)EMD(T)dT (1) 

We use the primary line (see Table[Q from every Fe ion in 
the observed spectra to get a set of integral equations (eq.Q, 
whose solution yields the EMD scaled by the unknown Fe 
abundance. For other elements, we use ratios of the He-like to 
H-like line fluxes instead of absolute fluxes, thus the element 
abundance A- cancels out. This adds another set of equations 
that constrain the shape of the EMD, and do not depend on the 
abundances: 



The X-ray spectra include as many as ten Fe ions but no 
more than two ions from other elements. In total, we get 14 
equations for the flare and quiescence observations, but dif- 
ferent equations, depending on which lines are visible (see 
Table [Q. We fit the line fluxes and flux ratios using the least 
squares best fit method to solve for the EMD, where the EMD 
is expressed by a parameteric non-negative function of T. This 
method yields the estimated shape of the EMD, independent 
of any assumptions for the abundances, and is scaled by the Fe 
abundance. The integration in eq.^and eq.[2]is cut-off at 8 keV 
beyond which the EMD is completely degenerate. This means 
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Table 1. Measured line fluxes used for EMD fitting and abundance calculations. 



Ion 


Wavelength 
A 


io- 4 


Flare 
photons citT 2 s 


Quiescence 
10~ 4 photons cm -2 s 


Ratio 
Flare/Quies. 


p e I6+ 


15.01 




6.61 ± 0.7 


4.85 ± 0.50 


1.36 ±0.20 


Fe 17+ 


14.20 




3.98 ± 0.44 


3.2 ± 0.33 


1.25 ±0.19 


Fe 18+ 


13.50 




4.08 ± 0.61 


3.1 ±0.32 


1.32 ±0.24 


Fe 19+ 


12.83 




4.98 ± 0.59 


3.08 ± 0.32 


1.62 ±0.26 


Fe 20+ 


12.28 




3.64 ± 0.53 


3.14 ±0.33 


1.16 ±0.21 


Fe 21+ 


11.77 




2.89 ± 0.52 


2.18 ±0.23 


1.33 ±0.28 


Fe 22+ 


11.00 




3.19 ±0.6 


1.40 ±0.15 


2.29 ± 0.49 


Fe 23+ 


10.64 




4.74 ± 0.63 


1.18 ±0.13 


4.02 ± 0.69 


Fe 24+ 


1.85 




2.38 + 0.3* 


0.136 ± 0.08 


17.6 ± 10.5 


Fe 25+ 


1.78 




0.51 +0.16* 






N 6+ 


24.78 


i 


5.66 ± 0.75** 


3.54 ± 0.46** 


1.88 ±0.32 


q6+ 


21.60 




2.37 ± 0.37 






o 7+ 


18.97 




21.8 ±2.2 


17.4 ± 1.8** 


1.25 ±0.18 


Ne 8+ 


13.45 




3.55 ± 0.75 


3.53 ±0.37 


1.01 ±0.24 


Ne 9+ 


12.13 




22.1 ±2.3 


16.0 ± 1.6 


1.39 ±0.20 


Mg 10+ 


9.17 




0.94 ± 0.27 


0.88 ±0.10 


1.07 ±0.33 


Mg u + 


8.42 




4.82 ± 0.66 


1.56 ±0.16 


3.08 ± 0.53 


Si I2+ 


6.65 




1.41 ±0.41 


1.14 ±0.12 


1.24 ±0.38 


Si I3+ 


6.18 




1.37 ±0.7 


1.12 ±0.12 


1.23 ± 0.64 


S 14+ 


5.05 






0.54 ± 0.07 




S 15+ 


4.73 






0.30 ± 0.05 




Ar 16+ 


3.96 






0.27 ± 0.04 




Ar 17+ 


3.73 






0.036 ± 0.03 




* Extracted from EPIC-pn 
** Used for abundances only 



NOTE: For He-like ions, the line flux includes the resonant line only. For H-like ions the 
flux includes both lines of the Isya doublet. For Fe the flux includes all the transitions 
of the ion within 0.03A of the given wavelength. Flux errors include 10% instrument 
calibration error. 



that some of the EM in the last bin could be attributed to even 
higher temperatures. 

The atomic dat a for the line powers ar e calculated using the 
HULLAC code dBar-Shalom et alll200ll) . In order to measure 
the line fluxes and solve for possible blending, we preform an 
ion-by-ion fitting to the spectra. The line powers for each ion 
are calculated at its maximum emissivity temperature and then 
passed through the instrument response. The observed spectra 
are fitted by a set of complete individual-ion spectra simulta- 
neously, resulting in an excellent fit that accounts for all the 
observed lines and b l ends. This process is similar to th e one 
used in iBehar. et al] fcOOll) and Brinkman et all feOOll) . The 
line fluxes used in the EMD fitting are listed in Table [2 The 
ionic abundance s ( f n ) for: Fe, Ar, S, Si , Mg are taken from lGu I 
J2003L whereas lMazzotta et al] fl998) is used for the other el- 
ements. 

Our goal is to compare the EMD of the flare and quies- 
cence states. It is important to note that the solution for the 
EMD is not unique as is the c ase with integral equations of 
this sort (Crai g & Brown! Il976l) . On scales much smaller than 
the width of the ions emissivity curves, or in temperature re- 
gions where there are no emissivity peaks of any ion, there is 
no way of constraining the EMD. Therefore, in order to be able 
to compare the EMD solutions, the confidence intervals of the 



solution are as important as the actual values. We choose to fit 
a staircase shaped function to allow for local confidence in- 
tervals estimates. The confidence intervals are calculated using 
the inverse x 1 distribution, meaning we search the parameter 
space for the x 1 contour that gives a deviation from the best fit 
that corresponds to the requested confidence level. 

The fitted EMD is plotted in Figure@]with 90% confidence 
intervals. The value of the EMD in each bin represents the aver- 
age EMD over the bin. Selecting the number of bins and their 
widths is not trivial. Since, as discussed above, we are inter- 
ested in meaningful confidence intervals, we cannot use nar- 
row bins, as this will result in excessive error bars. The line 
emissivity curves have considerable widths and some extend to 
temperatures much higher than their peak emissivity, resulting 
in strong negative correlations between the EM in neigbouring 
bins. The errors on the measured fluxes increase the uncertainty 
on the EM even further. Ultimately, if meaningful confidence 
intervals are to be obtained, the number of bins has to be kept 
small and their width optimized according to the constraints in 
each region. 

The important physical quantity is the integral of the EMD 
over a range of temperatures. The integrated EM from zero to 
kT is plotted in the bottom panel of Figure 0] with 90% confi- 
dence bars. The uncertainties caused by the strong correlation 
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Fig. 4. Top: EMD of the two observations. Error bars indicate 
90% confidence intervals, bottom: The integrated EM up to 
kT, with 90% confidence intervals. EMD is scaled according 
to Solar Fe abundance, taken to be: Fe/H = 4e-5. 

Table 2. The calculated abundances relative to Fe, obtained 
from the EMD analysis of the Chandra (quiescence) and XMM- 
Newton (flare). "Flare" abundances are in fact averaged "qui- 
escence+flare" components during the 2001 flare observed by 
XMM-Newton. Uncertainties are lcr. 





Flare 


Quiescence 


Ratio 


Element 


X/Fe 


X/Fe 


Flare/Quies. 


N 


15.1 ±2.7 


13.7 ±2.3 


1.1 ±0.3 


O 


30.3 ± 5 


42 ±6 


0.72 ±0.15 


Ne 


15.2 ± 2.3 


19.9 ± 2 


0.76 ±0.14 


Mg 


1.94 ± 0.33 


1.66 ±0.17 


1.17 ±0.23 


Si 


0.99 ± 0.4 


1.37 ±0.14 


0.72 ± 0.30 



between the EMD bins disappear with integration, resulting in 
much smaller error bars. 

In order to extract the X/Fe abundance ratios, we simply 
calculate the non-Fe line fluxes (eq. Q from the Fe-scaled 
EMD. The ratio between the measured and calculated flux 
gives the abundance value. The results for the elements ob- 
served by both instruments are summarized in Table [2] Errors 
for the abundances include both uncertainties in the EMD 
model and in the measured line fluxes. 

3. Discussion 

As seen in Figure |3j the relative line intensities in flare and in 
quiescence match very well for wavelengths longer than 12 A. 
For the most part, the residual is less than lcr. The fluctuations 
in the residual around the brightest lines at 19, at 13.5, and at 
12 A are caused by slight inaccuracies in the LSF of the RGS 
intensified by the strong steep peak. However, the residuals do 
average to zero. Around the 19 A line, that is by design of the 
scaling. Even at shorter wavelengths, the lines of the He-like 
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Mg 10+ (9.17 A) and of Si 12+ (6.65 A), which have maximum 
emissivity temperatures of 540 and 860 eV, respectively, do not 
show detectable variations. 

More evident changes in the spectrum start to appear be- 
low 12 A. The most conspicuous residuals are those due to 
lines of highly ionized iron: Fe 22+ and Fe 23+ between 10 and 
12 A, whose emissivities peak at 1.2 and 1.6 keV, respectively. 
Since the relative intensities of the lower temperature Fe lines 
do not change during the flare, this is clearly a temperature ef- 
fect rather than an abundance effect. Together with the dramatic 
change in the shape of the continuum, it indicates the appear- 
ance of a very hot component in the flare with little effect on 
the cooler plasma. This result is also seen from the EMD in 
Fig. El 

Another line that shows a significant change between flare 
and quiescence is the Lyo- line of Mg 11+ at 8.42 A, whose 
flux increases by a factor of 2.4 ± 0.4 during the flare, be- 
yond the overall 25% enhancement. This change too is due 
to the high-T flare as explained in the next paragraph. This 
comes to demonstrate that without clear identification of high- 
T lines, abundance effects can be easily confused with temper- 
ature changes, casting doubt on the reliability of some earlier 
results, which could not rely on line resolved spectra. The only 
other ion marginally in the RGS band that could probe tem- 
peratures above 1 keV is Si 13+ , whose lines suffer from a low 
signal to noise ratio. 

From the EMD comparison in Fig.[4l we conclude that with 
high certainty, a large, hot component above 3 keV has risen 
during the flare and did not exist during quiescence. In the 2 - 
4 keV region we get poor constraints on the flare EMD, due to 
the absence of lines with peak emissivity in this temperature re- 
gion from the RGS band. This intensifies the negative correla- 
tion between bins mentioned above and results in a large error. 
For the quiescence EMD, we have Si, S and Ar lines that allow 
for much better constraints in this region. At 2 keV and be- 
low, both EMDs seem to be similar, with the flare EMD being 
slightly higher. The quiescent LETG spectrum is missing lines 
from wavelengths longer than 20 A due to the aforementioned 
problems with CC mode. These lines would have allowed for 
much better constraints on the quiescence EMD below 1 keV. 

We conclude that the plasma in its quiescent phase 
(< 2 keV) was not affected by the flare, except its emission 
measure increased by a constant factor of ~1.25. We note that 
the large amount of high-T EM in the flare contributes some 
flux to lines of low-T H-like ions, since their emissivity ex- 
tends to higher temperatures, much beyond the peak emissiv- 
ity. This explains the added flux to the Mg 1 1+ (8.42 A) line. For 
the lines of the H-like ions 7+ and Ne 8+ , the flare component 
contributes ~ 30% of the total flux in the line. This is why our 
model requires that the O and Ne abundances decrease by that 
amount to maintain the relative quiescent flux level (Tables- 
No significant abundance variations beyond these are observed. 

This means that most likely, the 25% brightening of the 
corona is due to gradual variations in the quiescent corona over 
the 15 months that passed between observations. In 11 obser- 
vations by ROSAT over a similar time span, gradual changes 
of over 40% in the total count rate were registered dYietall 
119971 These variations may correspond to the activity cycle of 
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the system or changes in the active regions rather than to flares. 
Aud ard et alJ J200ll) have also reported that during a flare on 
HR 1099 the colder quiescent plasma was not affected. Both 
HR 1099 and <x Gem are constantly very bright X-ray sys- 
tems. Therefore, one expects to get a large emission contribu- 
tion from the corona outside the flaring region. 

Simple heating of the plasma does not change the total EM, 
therefore variations in the total EM indicate added material, 
changes in density or both. The total integrated EM up to 3 keV 
is slightly higher during the flare, although still consistent with 
the quiescence EM, within the 90% confidence intervals. The 
total EM up to the 8 keV cut-off during the flare is 5.4+1.9 
times that of the quiescence EM, which is very large in a bright 
system such as <x Gem, so it is unlikely that such a huge amount 
of plasma is added to the corona. The more likely interpretation 
is that the hot EM originates from plasma heated lower in the 
chromosphere, where the higher density would result in a large 
EM, even for a small amount of evaporated plasma. The in- 
crease in density was not detected here, but this could be due 
to the high charge states (H-like and bare) typical of the high-T 
flar e, for which no dens i ty diag nostics are available. 

ISanz-Forcada et aD J2£)02) analysed the 1998 flare on 
cr Gem, detected by EUVE. The flare timescale of approxi- 
mately 1 day, as well as the sharp rise and gradual decline are 
similar to those of the present flare. They report an increase in 
line flux of highly ionized Fe by a factor of 2-4 and a contin- 
uum increase by a factor of 6. Their EM analysis shows that an 
excess EM is prevalent throughout a broad temperature range. 
However, the absence of error estimates for the EMD does not 
allow for unambiguous conclusions. It is possible that the flare 
reported by Sanz-Forcada et al. was similar to the present one, 
only the EUV measurements were limited to probing temper- 
atures of up to 1.5 keV, which is below the lower temperature 
end of the present flare. In the bottom panel of Fig. 0] it can 
be seen that there is no appreciable effect on the integrated EM 
up to that temperature. The fact that no significa nt density in- 
crease was detected bv lSanz-Forcada et alJ (120021) may hint to 
the presence of much higher temperatures and ionization states 
in that flare, as were observed here. 

4. Conclusions 

There is little to no change in the chemical and thermal com- 
position under 2 keV between the flaring and quiescent states. 
This seems to suggest that the flare did not have a significant in- 
fluence on the coronal structure outside a limited flaring region. 
The effect of local flaring is similar to that which is observed 
on the sun. If the hot flaring comp onent we observe i s indeed, 
as suggested by the Neupert effect (Giidel et al., 2002), plasma 
evaporating from the chromosphere to the higher corona, then 
as it fills the coronal loop, the EM of the evaporated plasma will 
first increase in the loop as more material flows in, although 
it is already cooling, while in a later phase, it will cool back 
to chromospheric temperature and then condensate out of the 
corona, i.e., the EM will decrease again. During this process, 
most of the surrounding cooler plasma that we observe in the 
quiescent spectrum remains unaffected. As the flare decays, its 
contribution rapidly decreases and becomes insignificant to this 



quiescent emission. The fact that the large flare had little influ- 
ence on the cooler quiescent corona, suggests that the corona is 
stable against single large events and element enrichment pro- 
cesses are probably slow and continuous. 
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